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ABSTRACT 



Context. High-resolution soft X-ray spectroscopic observations of single hot white dwarfs are scarce. With the Chandra Low-Energy 
Transmission Grating, we have observed two white dwarfs, one is of spectral type DA (LB 1919) and the other is a non-DA of spectral 
type PG 1159 (PG 1520+525). The spectra of both stars are analyzed, together with an archival Chandra spectrum of another DA 
white dwarf (GD 246). 

Aims. The soft X-ray spectra of the two DA white dwarfs are investigated in order to study the effect of gravitational settling and 
radiative levitation of metals in their photospheres. LB 1919 is of interest because it has a significantly lower metallicity than DAs 
with otherwise similar atmospheric parameters. GD 246 is the only white dwarf known that shows identifiable individual iron lines in 
the soft X-ray range. For the PG 1 159 star, a precise effective temperature determination is performed in order to confine the position 
of the blue edge of the GW Vir instability region in the HRD. 

Methods. The Chandra spectra are analyzed with chemically homogeneous as well as stratified NLTE model atmospheres that assume 
equilibrium between gravitational settling and radiative acceleration of chemical elements. Archival EUV and UV spectra obtained 
with EWE, FUSE, and HST are utilized to support the analysis. 

Results. No metals could be identified in LB 1919. All observations are compatible with a pure hydrogen atmosphere. This is 
in stark contrast to the vast majority of hot DA white dwarfs that exhibit light and heavy metals and to the stratified models that 
predict significant metal abundances in the atmosphere. For GD 246 we find that neither stratified nor homogeneous models can fit 
the Chandra spectrum. The Chandra spectrum of PG 1520+525 constrains the effective temperature to r c fj = 150 000 ± 10 000 K. 
Therefore, this nonpulsating star together with the pulsating prototype of the GW Vir class (PG 1 159-035) defines the location of the 
blue edge of the GW Vir instability region. The result is in accordance with predictions from nonadiabatic stellar pulsation models. 
Such models are therefore reliable tools to investigate the interior structure of GW Vir variables. 

Conclusions. Our soft X-ray study reveals that the understanding of metal abundances in hot DA white dwarf atmospheres is still 
incomplete. On the other hand, model atmospheres of hydrogen-deficient PG 1159-type stars are reliable and reproduce well the 
observed spectra from soft X-ray to optical wavelengths. 
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1. Introduction Until recently, only five WDs were observed with Chandra 

LETG. H 1504+65 is an extremely hot (T eff = 200 000 K), pecu- 

An unexpectedly small number of white dwarfs (WDs) were liar PG 1 159 star ' obviously a naked C-O or O-Ne WD, whose 

detected in the ROSAT PSPC X-ray all sky survey (Fleming soft X " ra y spectrum is characterized by highly ionized O, Ne, 

et al. 1996). Most of them (161) are of spectral type DA (pure- and M 8 llnes ( Werner et aL 2004a )- ° n the other hand ' GD 246 

hydrogen optical spectra), with an additional three DOs (helium- 18 a DA ' and lts Chandra spectrum allowed the first unambigu- 

dominated), three DAOs (mixed H/He optical spectra), and eight ous identification of lines from highly ionized iron (Vennes & 

PG 1159 stars (He-C-O dominated). It was realized that the Du P ms 2002 )- The Chandra spectra of the two DAs HZ43 and 

atmospheric opacity of radiatively levitated metals effectively Slnus B are perfectly matched by pure-hydrogen atmospheres 

blocks the outward leakage of X-rays in the vast majority of DAs and "« used as soft-X-ray calibration targets (Pease et al. 2000; 

with r eff > 40 000 K (Barstow et al. 1997). Furthermore, the in- Vennes & Du P uls 2 002; Beuermann et al. 2006). A Chandra ob- 

terstellar medium proved denser than expected in many lines of servation of the hot DO KPD 0005+5 106 (r eff = 200 000 K) was 

sight (see, e.g., Welsh et al. 1999), thus affecting cooler WDs. used to P rove that the soft X-ray emission is of photosphenc and 

not coronal origin; however, the S/N ratio is insufficient to iden- 

High-resolution soft X-ray spectroscopy of hot WDs enables tify individual spectral lines (Drake & Werner 2005). 

the identification of chemical species in their photospheres that Few other WDs are bright enough in soft X-rays to ob- 

cannot be detected in other wavelength ranges. However, such tain useful Chandra spectra. We have observed two of them 

observations are rather scarce and only became feasible with (LB 1919 and PG 1520+525) and including archival data of 

the advent of the Chandra observatory with its Low-Energy GD246 present a spectral analysis in this paper. The DA 

Transmission Grating (LETG). LB 1919 is of interest because previous Extreme Ultraviolet 
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Explorer (EUVE) spectroscopy showed an unexpectedly high 
soft X-ray flux that reveals an unexplained low metal abundance 
(Wolff et al. 1998), in contrast to expectations from radiative - 
levitation theory. We performed X-ray spectroscopy in order to 
identify individual elements that could hint at the origin of the 
metal deficiency. 

We observed PG 1520+525 in order to constrain its effective 
temperature. The motivation lies in the fact that this nonpulsating 
PG 1159 star, together with the pulsating prototype PG 1 1 59— 
035, defines the blue edge of the GW Vir instability region in the 
Hertzsprung-Russell diagram (HRD), see e.g. Jahn et al. (2007). 
The homogeneous atmosphere of PG 1520+525 provides a con- 
venient comparison for the investigated DA white dwarfs be- 
cause it is not subject to the uncertainties in the treatment of ra- 
diative levitation and stratification of atmospheric composition. 
Thus, PG 1520+525 provides a test case for the model atmo- 
spheres employed in this study. A failure to achieve a consis- 
tent fit to the UV/optical and X-ray spectra of the homogeneous 
PG1159 photosphere would indicate serious model deficiencies 
in addition to problems that would be encountered with radiative 
levitation physics in the case of hot DA WDs. 

This paper is organized as follows. We first describe in more 
detail the motivation of our soft X-ray analyses (Sect. 2). Then 
we specify our model atmosphere calculations and the atomic 
data used (Sect. 3). In Sect. 4, we report on the observations uti- 
lized in our analysis. In Sections 5 to 7 we delineate the analysis 
procedure of our three program stars one at a time, and in Sect. 8 
we conclude. 



2. Soft X-ray emission from white dwarfs 

A concise presentation of past soft X-ray observations of WDs 
and their interpretation can be found, e.g., in Werner (2008). 



2.1. DA white dwarfs 

DA WDs possess almost pure hydrogen atmospheres. The high 
gravity in these objects results in a chemical stratification with 
the lightest element, hydrogen, floating atop. In hot DAs (T e s > 
20 000 K), hydrogen is almost completely ionized and the opac- 
ity in the atmosphere is strongly reduced. Soft X-ray radiation 
can emerge from deep, hot, photospheric layers. However, this 
radiation may be blocked by the opacity of heavy elements that 
can be kept in the atmospheres by radiative levitation (Chayer 
et al. 1995a,b). The increasing efficiency of this mechanism 
with T e g results in only a few X-ray detected DAs with T e g 
> 60000 K. 

Since the spectral resolution of EUVE is too low to identify 
lines of individual species, the metal abundances of the hitherto 
investigated DAs were determined relative to the well-studied 
G 191-B2B (Wolff et al. 1998) using chemically homogeneous 
model atmospheres. The scaling of the metallicity in these mod- 
els relative to that of G 191-B2B (the so-called metallicity in- 
dex) results in satisfying spectral energy distribution (SED) fits 
for many, but not all, DAs. 

To make further progress, non-LTE model atmospheres were 
developed that calculate the abundances of the elements at each 
atmospheric depth in a self-consistent way, assuming equilib- 
rium between gravitational downward pull and radiative upward 
acceleration (Dreizler 1999). The EUVE SEDs of most DAs can 
be reproduced well by these models (Schuh et al. 2002). Some 
stars show a higher metallicity than predicted. This can be ex- 
plained, for instance, by ongoing accretion of circumstellar or 
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Fig. 1. Top panels: Element abundances over depth following from 
two different treatments of radiative acceleration by photoionization 
(Jeff = 54 000 K, logg = 8.2). Left: the whole photon momentum is 
transferred to the ion. Right: only a fraction is transferred, according to 
Eq. 3. Bottom panel: difference. The peaks correspond to a shift of the 
steep regions of the abundance pattern by one depth point in the model 
atmosphere. 



interstellar matter. Other stars have a lower metallicity than pre- 
dicted by these models, and the reason for this is not known. 
A drastic example is HZ 43 A (r e jf = 51 000 K, Beuermann et al. 
2006) whose SED is perfectly matched by a pure H model atmo- 
sphere. Another low-metallicity DA is the even hotter LB 1919 
(T e « = 56 000 K, this paper), which is subjected to analysis here. 

The exact knowledge of element abundances in the atmo- 
spheres of the metal-poor DAs might give hints as to an expla- 
nation for their metal deficiency. Many spectral lines of the high 
ionization stages of heavy elements like iron and nickel lie in 
the soft X-ray wavelength range. Spectroscopy with Chandra en- 
ables us to investigate these objects more closely. Our work rep- 
resents the first detailed analysis of high-resolution soft X-ray 
spectra of DAs. Besides LB 1919, we study the more metal-rich 
DA GD 246 as a reference object. 

2.2. PG 11 59 stars 

As mentioned in the introduction, only very few non-DAs (DO 
and PG1159 stars) were detected in the soft X-ray band be- 
cause of the additional opacity of helium and enriched metals. 
PG 1 159 stars are hot, H-deficient (pre-) WDs with atmospheres 
mainly composed of He, C, and O. It is thought that their H- 
envelope was consumed during a late He-shell flash (Werner & 
Herwig 2006, and references therein). Among the PG 1 159 stars, 
PG 1520+525 (r eff = 150 000 K) is one of the brightest soft X- 
ray sources and hence an interesting target to study the spec- 
tral characteristics of this class in that wavelength region. The 
surface abundances of these hot stars are not affected by gravi- 
tational and radiative acceleration because of a weak radiation- 
driven wind (Unglaub & Bues 1997). Consequently, their atmo- 
spheres and envelopes can be assumed to be chemically homo- 
geneous. 

Our primary aim for the Chandra observation of the non- 
variable PG 1520+525 is to constrain its effective temperature. 
Then we can compare its position in the log r e (j-log g diagram 
with that of the prototype of the GW Vir pulsators, PG 1159- 
035. The blue edge of the instability region is confined by these 
stars (Werner et al. 1996). Based on nonadiabatic computations, 
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Fig. 2. Effect of modified treatments of acceleration by line and photoionization transitions on occupation numbers of Fe v and Fe vi superlevels 
(T tf! = 56000 K, logg = 7.9). The six panels show: 1. Original treatment of bb and bf transitions. 2. New treatment of bf transitions. 3. Difference 
between the first two. 4. New procedure for bb transitions. 5. bf and bb refinements are implemented. 6. Difference between the last two. The 
legend in panel 1 refers to the superlevel designations. 



the most advanced pulsation models for GW Vir stars, are those 
presented by Quirion et al. (2007). Accordingly, the exact lo- 
cation of the blue edge depends on the envelope composition, 
primarily on the C and O abundance and, to a smaller extent, on 
the metallicity. Therefore, the concept of a blue edge is neces- 
sarily "fuzzy". It was shown by spectroscopic analyses that both 
stars have (within error limits) the same abundances of C, O, and 
Fe (see below) and thus define the blue edge for that particular 
composition. We can therefore prove or disprove the prediction 
of corresponding pulsation models as to the location of the edge. 
This represents a strong test for these models and their potential 
to derive the interior structure of PG 1159 stars by asteroseis- 
mology methods. 



3. Model atmospheres and atomic data 

For the spectral analysis, chemically homogeneous NLTE 
model atmospheres were computed with the Tubingen Model 
Atmosphere Package (TMAP 1 , Werner et al. 2003), and chem- 
ically stratified models with a special variant, the diffu- 
sion/levitation NGRT code (Dreizler 1999). In comparison to 
earlier work with NGRT, refinements to the code were applied 
for more realistic physics. Moreover, new model atoms for Fe, 
Ni, and Ge were constructed. These improvements are described 
in the following. 



3.1. Improvements to the diffusion/levitation code NGRT 

In the NGRT code, the element abundance pattern at any depth 
point in the atmosphere is calculated self-consistently by assum- 
ing equilibrium between gravitational and radiative forces. As a 
result, the atmosphere is no longer chemically homogeneous but 
vertically stratified. 



3.1.1. Bound-free transitions 

The primary assumption for the calculation of the radiative force 
of a bound-free (bf) transition in NGRT is that the photon mo- 
mentum is transferred completely to the ion, keeping the elec- 
tron out of consideration. As a consequence, the calculated ra- 
diative acceleration of the remaining ion might be overestimated. 
To check this assumption quantitatively, a correction according 
to Gonzalez et al. (1995) was implemented and tested. 

The radiative force on a particle of element A, with ionic 
charge i for a bound-bound (bb) or bf transition j, caused by 
photons in the frequency range (v, v + dv) is 



F'jdv - '^-dv = crJv)—dv, 
dt c 



(1) 



with the net momentum dp v , transported by the radiation flux 
T v . As mentioned above, in the case of a bf transition, the mo- 
mentum of the photon causing an ionization of ion A 1 is not 
completely transferred to ion A 1+l . A part of it is taken away 
by the ejected electron. Thus, a correction factor f lon has to be 
introduced, describing the remaining fraction of the momentum 
transferred to the ion A l+l . 

The correction factor can be expressed in terms of the fre- 
quency v of the photon, the threshold frequency vo necessary to 
eject the electron and a factor ai. 



4 (v - v ) 

/ion = I - a\. 

3 v 



(2) 



For testing purposes, the quantum calculation of Sommerfeld 
(1939) was taken and applied to the code. With his value of 
a\ =6/5, /i on becomes 



/ion 1 



% hv- X 
5 hv ' 



(3) 



http : // astro . uni - tuebingen . de/ - TMAP 



with the ionization energy threshold x- An interesting result of 
Eq. 3 was already stated by Michaud (1970). For certain fre- 
quencies the electron can be ejected with more momentum than 
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Fig. 3. Grotrian diagram of Geiv. The resonance doublet is labeled. 

brought in by the photon (/ e i ec = 1 —/ion > 1). As a consequence, 
the ion will be pushed back in the atmosphere by the photoion- 
ization. Thus, including a detailed bf absorption calculation can 
in principle lead to either higher or lower outward directed ra- 
diative forces on an ion. 

The comparison between abundances at each depth of the 
atmospheres calculated with old and new treatments of bf tran- 
sitions reveals, however, only minor and practically unimportant 
deviations (Fig. 1); there is no effect on the SED, in particular on 
the spectral lines, of the model. 

3.1.2. Bound-bound transitions 

In the original NGRT code as used by Schuh et al. (2002), the re- 
distribution of transferred photon momentum over the ionization 
stages A 1 and A 1+1 by a bb transition in A 1 was treated such that 
the momentum was completely transferred to the next higher ion 
A 1+1 . This treatment is based on the assumption that an ionization 
following the radiative excitation of the ion takes place before 
the ion is deexcited by an inelastic electron collision. Especially 
in very hot atmospheres, this need not be true. For a realistic 
treatment, the probabilities for ionization and collisional deexci- 
tation would have to be evaluated. In a statistical sense, fractions 
of the momentum are transferred to the next higher ion A 1+1 and 
to the originally excited ion A 1 . 

As a rough approach to test the maximum possible effect, 
the calculation of the bb transitions was modified such that the 
entire photon momentum is transferred to the lower ion A 1 . For 
example, we show the effect of the modified bb (and bf) proce- 
dures on the population numbers of iron superlevels (introduced 
below) in Fig. 2. While the modified bf transitions leave the oc- 
cupation numbers essentially unaltered, the modification of the 
bb transitions results in an obvious change, although only in the 
outer photospheric layers. The effect on the emergent flux is be- 
low 10% and thus considered insignificant, even under the ex- 
treme assumption of the fully inverted momentum distribution. 
Therefore, a more complex treatment is unnecessary for our pur- 
poses. 

3.2. Model atoms 

For our model atmosphere calculations, we used two different 
sets of model atoms for the analysis of the PG 1 159 star on the 
one hand and for the two DAs on the other hand. The first set 
is identical to that described in detail by Werner et al. (2004a) 
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Fig. 4. GelV line profiles from LTE (dashed) and NLTE (solid) H+Ge 
model atmospheres (r cff = 57 000 K, logg = 7.6, Ge/H = 2.5 • 10~ 9 ). 
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Fig. 5. Comparison of H+Fe model spectra calculated with old and new 
Kurucz line lists (top left and right panels, respectively; T e g = 57 000 K, 
logg = 7.90, Fe/H = 1.25-10~ 6 ). Bottom panels: same models after con- 
volution with the Chandra LETG instrument response (bin size 0.1 A). 



in their Chandra analysis of H 1504+65. It comprises He, C, O, 
Ne, and Mg. The second set incorporates H, He, C, O, Si, P, S, 
and Ge, being represented by detailed model atoms taken from 
the Tubingen Model Atom Database 2 TMAD. It is augmented by 
Fe and Ni which are treated in a statistical way using superlevels 
and superlines (Rauch & Deetjen 2003). Table 1 summarizes the 
set of data used for the DA analyses. 

The Ge model atom is new and presented in more detail 
below. We also use new Fe and Ni model atoms that comprise 
vastly extended line lists (Kurucz 2009) compared to all previ- 
ous analyses. Their characteristics are also described in detail 
below. 

3.2.1 . Germanium 

Germanium (Ge iv) was first identified in three hot DAs, among 
them GD 246, by Vennes et al. (2005). For GD 246, we perform 
the first NLTE analysis of Ge in a stellar atmosphere and in- 
vestigate its radiative levitation properties. The model atom was 
constructed using level energies from the NIST 3 database. Due 
to the lack of data, all oscillator strengths fij are approximated 
by adopting values from the isoelectronic C iv, with the excep- 
tion of the two observed Ge iv AA 1 189.07, 1229.84 A lines, for 
which we use the /ij from Chayer et al. (2005). Very recently, 
new /ij for Ge iv lines were published (Nath Dutta & Majumder 

2 http : //astro . uni - tuebingen . de/ ~ TMAD 

3 http://physics.nist.gov/ 
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Table 1. Statistics of model atoms used in the calculations for 
LB 1919 and GD 246. Numbers in brackets denote levels and lines 
combined into superlevels and superlines. 



Elem. 


Ion 


NLTE levels 




lines 




H 


I 


10 




45 






II 


1 








C 


II 


16 




37 






III 


58 




329 






IV 


54 




295 






V 


1 








O 


II 


16 




26 






III 


21 




38 






IV 


18 




39 






V 


30 




95 






VI 


54 




291 






VII 


1 








Si 


II a 


20 




48 






III 


17 




28 






IV 


16 




44 






V 


15 




20 






VI 


1 








P 


III 


3 











IV 


15 




9 






V 


18 




12 






VI 


1 




_ 




S 


III 11 


3 











IV 


17 




32 






V 


19 




32 






VI 


18 




48 






VII 


1 




_ 




Ge b 


III 


14 











IV 


8 




8 






V 


9 











VI 


1 




_ 




Fe 


III 


7 


(2467) 


27 


(537689) 




IV 


7 


(6389) 


27 


(3102371) 




V 


7 


(6728) 


25 


(3266247) 




VI 


8 


(5464) 


33 


(991935) 




VII 


9 


(2690) 


41 


(200455) 




VIII 


8 


(694) 


33 


(19587) 




IX 


1 








Ni 


III 


5 


(3739) 


14 


(1033920) 




IV 


7 


(5939) 


27 


(2512561) 




V 


7 


(6165) 


25 


(2766664) 




VI 


7 


(10985) 


27 


(7408657) 




VII 


8 


(11866) 


33 


(4195381) 




VIII 


7 


(6469) 


27 


(1473122) 




IX 


1 









LB 1 9 1 9 models only b GD 246 models only 



201 1). A single test calculation was performed in which we use 
the new data for all line transitions in the model atom. The effect 
on the observed two lines is insignificant. Photoionization rates 
were computed with hydrogen-like cross sections. Electron col- 
lisional excitation and ionization rates were evaluated with usual 
approximation formulae following van Regemorter (1962) and 
Seaton (1962), respectively. 

In order to check for the importance of NLTE effects, we per- 
formed another test calculation in which we enforced LTE popu- 
lations for Ge levels by artificially increasing electron collisional 
rates. The NLTE line profiles are significantly shallower than the 
LTE profiles (Fig. 4). As a consequence, the derived NLTE abun- 
dances will be roughly a factor of two higher compared to LTE 
analyses. 




150 155 160 165 170 150 155 160 165 170 
X [A] X [A] 



Fig. 6. Like Fig. 5, but here for nickel (Ni/H = 7.9 ■ 10~ 8 ). 



Table 2. Observation log of spectra analyzed in our study. 



Star 


Instrument (Id.) 


Observation 
start date (UT) 
mm-dd-yyyy 


^exp 

(ks) 


LB 1919 


EUVE 


04-16-1994 


166 




FUSE LWRS 


05-03-2001 


33 




IUE SWP 52808 


11-15-1994 


4.5 




Chandra LETG/HRC-S 


02-02-2006 


111 


GD246 


EUVE 


08-08-1994 


14 




FUSE LWRS 


11-12-2000 


1.5 




FUSE MDRS 


07-14-2001 


24 




HST STIS/E140H 


11-20-1998 


2.4 




IUE SWP 40467 


12-27-1990 


0.33 




IUE LWP 19485 


12-27-1990 


0.84 




Chandra LETG/HRC-S 


01-14-2000 


40 


PG 1520+525 


Chandra LETG/HRC-S 


04-04-2006 


142 



3.2.2. Iron and nickel 

The Iron Opacity Interface (IrOnIc) constructs model atoms 
for iron-group elements as input for the stellar-atmosphere code 
(Rauch & Deetjen 2003). During the course of this work, new Fe 
and Ni line data were provided by Kurucz (2009) 4 . They include 
a much larger number of lines (about 20 times more), especially 
in the X-ray range. Figures 5 and 6 demonstrate the effect of 
these additional lines used in the atmosphere calculations. The 
new line lists were not incorporated in the bulk of the model cal- 
culations, however, the best-fitting models presented here were 
recalculated with the newest Kurucz atomic data sets available. 

4. Observations, reddening, interstellar H and He 
column densities 

We performed Chandra spectroscopy of LB 1919 and 
PG 1520+525 with the LETG/HRC-S. For GD246, we 
used the Chandra observation investigated by Vennes & Dupuis 
(2002). Other observations used in our work were retrieved 
from the MAST 5 archive. All observations are listed in Table 2. 

In order to determine interstellar reddening, FUSE, IUE, and 
HST-STIS spectra (the latter marked with "S" in Fig. 7) as well 

4 http : //kurucz . harvard . edu/atoms . html 

5 http://archive.stsci.edu/ 
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Fig. 7. Determination of interstellar reddening of LB 1919 and GD246 
from the observed flux distribution compared to models (thick red 
lines). For clarity, all spectra are convolved with 2 A (FWHM) 
Gaussians. 



as optical and IR magnitudes from the SIMBAD and NOMAD 
databases were investigated. Fluxes of our final models were 
scaled to the reddest available photometric magnitude. The red- 
dening law of Fitzpatrick (1999) was used with R v = 3.1. We de- 
termined £ B -v = 0.030 + 0.005 and 0.001 1+0.0001 for LB 1919 
and GD 246, respectively. These quantities will be used below 
for distance determination. Accounting for absolute-flux uncer- 
tainties in the UV, more realistic errors are of the order +0.05. 
However, the error in the distance determination is dominated by 
the uncertainties in T e ff and logg, and we neglect the errors for 
reddening. For PG 1520+525, Dreizler & Heber (1998) derived 
£ B -v = 0.0. 

Column densities for interstellar neutral hydrogen as well 
as neutral and ionized helium affect the analysis of EUVE and 
Chandra spectra. 

LB 1919: Based on the values found in the EUVE analy- 
sis by Landenberger-Schuh (2005) (N(Hi) = 1.6 • 10 19 ctrT 2 , 
N(Hei) = 0.04 N(Hi), NfHen) = 0.04 N(Hi)), the interstel- 
lar H and He column densities were varied in our fitting pro- 
cedure. Our values for two of the best-fitting models are de- 
picted in Fig. 12. For our present fit to the Chandra spectrum, 
we used N(Hi) = 1.6 ■ 10 19 cirT 2 , N(Hei) = 0.02 ■ N(Hi), 
N(Hen) =0.04-N(Hi). 

GD246: The Hi column density of 1.288 ■ 10 I9 ctrT 2 was 
taken from Oliveira et al. (2003) and kept fixed. Our adopted 
values for the column densities of Hei (0.12 N(Hi)) and Hen 
(0.02 N(H i)) are close to the ones determined by Vennes et al. 
(1993) (N(Hei) = 1.05 - 1.25 • 10 18 ctrT 2 , N(Hen) = 3.4 - 4.0 • 
10 17 cnT 2 ). 

Since these numbers depend on their assumed photospheric 
model, we adjusted the values slightly to better fit our mod- 
els. For our Chandra analysis presented here we chose 0.05 and 
0.026 N(Hi), respectively, which is close to the values used by 
Landenberger-Schuh (2005). 

PG 1520+525: Werner et al. (1996) used N(Hi) = 1.5 ■ 
10 20 crrr 2 for their fit of a T eff = 150 000 K model to the SED ob- 
served by EUVE. Dreizler & Heber (1998) determined N(H i) = 
2.5 TO 20 cirT 2 from the Lyo- line profile. For the fit to the Chandra 
spectrum, the H i column density was treated as a free parameter. 
For our final model fits in this paper, we found 1.0 • 10 20 cm~ 2 . 
The variation of the He column densities had no effect on the 



fit quality, and were kept fixed at N(He i) = 0.05 N(H i) and 
N(Hen) = 0.026 • N(Hi). 

5. LB1919 

5.1. Previous investigations 

LB 1919 (WD 1056+516) was analyzed by Vennes et al. (1997) 
using Balmer-line spectroscopy and pure-hydrogen model atmo- 
spheres. They obtained T e s = 68 640 K and logg = 8.08 (cgs 
units). Finley et al. (1997) excluded the star from their analysis. 
Their spectrum exhibited flat-bottomed Balmer-line profiles, and 
they suggested orbital or rotational velocities of ~ lOOOkm/s as 
a possible interpretation, although they may be artifacts (Koester 
priv. comm.). More recently, Lajoie & Bergeron (2007) deter- 
mined r e ff = 67 022 K and logg = 7.99, and Gianninas et al. 
(201 1) found r eff = 68 5 10 K and log g = 7.94 by pure-hydrogen 
model fits to the Balmer lines. 

Wolff et al. (1998) investigated a sample of 20 DAs from 
the EUVE archive. Introducing the above-mentioned metallicity 
index (Sect. 2.1), they were able to reproduce the spectra of most 
of them. For LB 1919, they determined a metallicity at least ten 
times smaller compared to that of G 191-B2B, adopting and 
logg values from Vennes et al. (1997). 

Based on these results, Landenberger-Schuh (2005) per- 
formed an analysis of EUVE spectra of DAs using new chemi- 
cally stratified NLTE model atmospheres with equilibrium abun- 
dances. In the case of LB 1919, the use of these new models re- 
sulted, however, in a worse fit compared to the analysis of Wolff 
et al. (1998) with homogeneous models because of too much 
levitation of metals in the diffusion models. The reason for this 
overprediction remained unexplained. 

5.2. FUSE 

A grid of pure-hydrogen model atmospheres was used to fit 
the Lyman lines in the FUSE spectral range. We found 7^ = 
56 000 + 2000 K and logg = 7.9 + 0.3. Somewhat hotter mod- 
els with higher gravity (e.g., T e g = 62 000 K, \ogg = 8.2) 
fit almost equally well, but are excluded by the Chandra spec- 
trum (Sect. 5.4). The temperature is significantly lower than 
the Balmer-line results reported above. It is remarkable that 
Lajoie & Bergeron (2007) arrived at a similarly low value of 
57701 + 18400K from the overall UV (ifJE)/optical flux dis- 
tribution, and at an even lower value of 39 045±18 000K by 
fitting the Lya line in the IUE spectrum. Figure 8 shows the 
comparison between a hotter model and a better-fitting, cooler 
model together with the FUSE Lyman lines. Adding metals to 
the 56 000 K model with abundances at the upper limit derived 
below has no significant influence on the Lyman-line profiles. 

Upper limits to metal abundances were determined with ho- 
mogeneous model atmospheres. Their accuracy is estimated to 
be about 0.3 dex following from the error ranges in 7^ and log g, 
and they are given in number ratios relative to H. We also com- 
pare line profiles from diffusion models with observations. 

Carbon and Oxygen From the absence of Cm Al 175 A, we 
derive an upper limit of C/H = 1.0 • 10~ 7 . The stratified models 
overpredict the line strengths (Fig. 9). Note that all computed 
line profiles have zero rotation velocity. 

Except for the Ovi AA 1031.91, 1037.61 A resonance dou- 
blet, no O line is identified. Because of a possible ISM contri- 
bution, only an upper limit of O/H = 1 • 10~ 6 can be derived. 
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Fig. 8. Fit to the Lyman lines of LB 1919 (thin line). Overplotted are 
pure H models with T^g = 56 000 K, log g = 7.9 (thick line), and with 
T cff = 69 000 K, logg = 7.9 (dashed line). 
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Fig. 9. Homogeneous (thick line) and stratified (dashed line) model 
spectra near the Cm multiplet at A 1175 A in LB 1919. The homoge- 
neous model has C/H = 1.0 • 10~ 7 . 
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Fig. 10. Homogeneous (thick line) and stratified (dashed line) model 
profiles of P v and Si iv compared to LB 1919. The homogeneous model 
has Si/H = 5.0 • 10" 9 and P/H = 1.0 ■ lO" 10 . 



The stratified models fit the observation best for logg = 8.2. A 
smaller gravity results in too strong features. 



Silicon, Phosphorus and Sulfur From the absence of Siiv 
AA 1066.63, 1 122.48, 1 128.33 A, P v AA 1 1 17.98, 1 128.01 A, and 
Siv 1062.66, 1072.96, 1073.508 A, upper limits of Si/H = 
5.0 • 10~ 9 , P/H = 1.0 • 10~ 10 and S/H = 1.75 • 10~ 8 are de- 
rived (Figs. 10, 1 1). In all cases, the diffusion models predict too 
strong line profiles. 

Iron and Nickel No lines are identified corresponding to upper 
limits of Fe/H < 1.0 • 10~ 6 and Ni/H < 1.0 • 10~ 6 . At T ef[ = 
56 000 K, the stratified models predict no detectable lines only if 
the surface gravity is logg = 8.2 or higher. 



2.25 
2.00 

x 1 " 75 
| 1-50 
j? 1.25 
1.00 
0.75 



T, B = 56 000 K 




1063.5 1066.5 1067.5 1072.5 
MA] 



Fig. 11. Homogeneous (thick line) and stratified (dashed line) model 
profiles of S iv and Si iv compared to LB 1919. The homogeneous model 
has S/H = 1.75 • 10~ 8 and Si/H = 5.0 • 1Q- 9 . 
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Fig. 12. Best-fit homogeneous (solid) and stratified (dashed) models for 
the EUVE spectrum of LB 1919 (thin line). The homogeneous model 
has metal abundances set to the upper limits given in Table 4; neither 
model has oxygen. 



5.3. EUVE 

The stratified models calculated by Landenberger-Schuh (2005) 
had a flux too low to match the EUVE observation. The failure 
to derive a good fit resulted from the excessive opacity caused 
by the large number of elements incorporated in the model at- 
mospheres. We tried to reproduce the EUVE observation with 
stratified models that include only those metals listed in Table 4, 
i.e., metals that have been detected in the FUSE spectra of com- 
parable WDs, like GD 246. Our results show that the strongest 
influence on the EUV SED comes from the opacity of oxygen 
and that a satisfying fit can only be derived for a model atmo- 
sphere without this species and at high gravity (logg = 8.5). We 
also computed homogeneous models with the same elements and 
abundances set to the upper limits derived from the FUSE analy- 
sis (again without oxygen) and also obtained a good fit (Fig. 12). 
The same result is obtained from a pure H model. 



5.4. Chandra 

No Fe or Ni lines are detected in the Chandra spectrum, giving 
upper abundance limits of Fe/H and Ni/H < 1.25 ■ 10~ 7 derived 
from homogeneous Fe+H and Ni+H models. These limits are 
almost one order of magnitude more stringent than the FUSE- 
derived results. The observed overall SED can be reproduced 
successfully with both homogeneous and stratified models. The 
interstellar H and He column densities were kept fixed to the 
El/VE-derived values. For the stratified models, however, there 
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Fig. 13. Best-fit models (thick lines) for the Chandra spectrum of 
LB 1919 (thin lines). Top: homogeneous model, bottom: diffusion 
model. Both models are computed with H, C, O, Si, P, and S. The ho- 
mogeneous model has abundances equal to the upper limits given in 
Table 4. 



is still enough Fe and Ni levitated that lines should be detectable 
in the Chandra spectrum, even at a higher gravity of log g = 8.5. 

For further analysis, homogeneous and stratified models 
were calculated that incorporated the same metals as the EUVE 
analyses in the preceding section. Individual spectral features in 
the observation are not recognizable, which serves as an addi- 
tional constraint for the upper limits of the metal abundances de- 
termined in the FUSE analysis. Both kinds of models match the 
Chandra spectrum well. The models displayed in Fig. 13 include 
oxygen, the homogeneous model with an abundance according 
to the upper limit given in Table 4. The vast majority of the weak 
lines in these models stems from O. In essence, a pure H model 
gives an equally good fit. 

5.5. Mass and distance 

To determine the mass from T e s and log g, we use DA evolution- 
ary tracks by Althaus & Benvenuto (1998). We chose the tracks 
with metallicity z = 10 3 and hydrogen-envelope fractional mass 
of Mh/M* = lfT 4 and derive M = 0.66 !^M 3 . Taking the flux 
calibration of Heber et al. (1984), we compute the distance, us- 
ing the extinction-corrected visual magnitude that is determined 
from the observed magnitude V= 16.41 to my = V - 2,11 '4- 
E(B-V) with E(B-V) = 0.030. With r eff = 56 000 + 2000 K, 
logg = 7.9 + 0.3, and H v = 7.33 • 10~ 4 erg cirT 2 sec" 1 Hz" 1 
(model Eddington flux at 5454 A) we find (see, e.g., Rauch et al. 
1994) 



d = 7.1 1 ■ 10 4 A /// v M10 a4mi 'o- 1 °ss = 321 1 



101 DC 

-137 P L ' 



(4) 



5.6. Summary on LB 191 9 



Our fit to the Lyman lines and the Chandra spectrum gave 
r eff = 56 000 + 2 000K, logg = 7.9 ± 0.3. The temperature 
is significantly lower (w 100007T) than that derived from pre- 
vious Balmer-line fits. This result is surprising because Barstow 
et al. (2003a) had established that Lyman-line analyses for DAs 
exceeding T e ff « 50 000 K yield systematically higher tempera- 
tures than Balmer-line analyses, i.e., the opposite that we found 
for LB 1919. In their sample, however, there is also one excep- 




Fig. 14. Results from previous analyses of GD246 (Table 3). Results 
obtained by LTE analyses (open circles) are compared to NLTE anal- 
yses (filled circles) and the analysis with stratified NLTE model atmo- 
spheres (filled square). The triangles mark the two results of this paper, 
obtained by fitting the Lyman lines (full symbol) and the Chandra SED 
(open symbol). 



tion to this general trend. Very similar to LB 1919, PG 1342+444 
has r e ff = 66 750 K and 54 308 K derived from Balmer and 
Lyman lines, respectively (see also Barstow et al. 2002). 

Metal lines are not detected in the FUSE and Chandra spec- 
tra. We derived upper abundance limits with homogeneous mod- 
els. The respective models fit the SEDs observed by Chandra and 
EUVE. Stratified models, on the other hand, overpredict metal 
abundances. Reducing the abundances in these models to bring 
them into agreement with observations requires an excessively 
high gravity (logg > 8.5), which is not compatible with the 
Balmer and Lyman lines. In essence, this confirms the result ob- 
tained from the EUVE analysis by Landenberger-Schuh (2005). 
All observations are compatible with a pure H atmosphere. 



6. GD246 

6.1. Previous investigations 

Earlier determinations of T e ff and logg of GD246 
(WD 2309+105) are summarized in Table 3 and Fig. 14. 
The most recent study by Gianninas et al. (2011) is probably 
the most reliable NLTE Balmer-line analysis because it utilizes 
new Stark broadening data (Tremblay & Bergeron 2009). (We 
use the same data in our analysis.) All analyses are based on 
Balmer- and/or Lyman-line profiles with the exception of two, 
which are based on the EUVE SED. One of them (Schuh 
et al. 2002) uses stratified NLTE models, and it results in an 
exceptionally large gravity. 

A small number of lines of trace metals could be identified in 
the FUSE and HST spectra of GD 246. Beside elements known 
to appear in WDs like Si, C, and P, Vennes et al. (2005) identified 
for the first time lines of germanium (Ge iv) in a HST observa- 
tion and determined log (Ge/H) = -8.6 ± 0.2. 

6.2. FUSE and HST 

Pure hydrogen models were used to fit the Lyman lines. We 
found r eff = 57 000 ± 2000 K and logg = 7.6 ± 0.3 (Fig. 15). 
Also plotted are the line profiles of one of our models that has 
parameters that are essentially equal to those found by Barstow 
et al. (2003a) from their Lyman-line analysis (T e s = 52 000 K, 
logg = 7.9). It can be seen that the differences between the two 
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Table 3. T cff and logg of GD 246 as determined by previous analyses and the present work (see also Fig. 14). 



No. 


T a TK1 
1 eff L^J 


lnc q I cm 1 v^l 


Method and Model Atmospheres 


Authors 


1 


60100 


7.72 


T TP 
Li L 


B aimer 


rl 


Vennes et al. (1997) 


2 


58 700 


7.81 


T TP 
LI L 


B aimer 


rl 


Finley et al. (1997) 


3 


59 000 


7.80 


T TP 
LI L 


JZU VIZ 


H-H metals 


Wolff etal. (1998) 


4 


53 100 


7.85 


TsJT TP 


B aimer 


1 1 
1 1 


Napiwotzki et al. (1999) 


5 


54000 


7.80 


LTE 


Lyman 


H 


Wolff etal. (2001) 


6 


56 000 


8.20 


NLTE 


EUVE 


H+metals stratified 


Schuh et al. (2002) 


1 




1 CO 


NLTE 


Lyman 


H+metals 


rsaistow et al. {zuujsi) 


8 


51300 


7.91 


NLTE 


B aimer 


H+metals 


Barstow et al. (2003a) 


9 


54400 


7.90 


NLTE 


B aimer 


H 


Liebert et al. (2005) 


10 


57 007 


7.82 


LTE 


B aimer 


H 


Koester et al. (2009) 


11 


54150 


8.00 


NLTE 


B aimer 


H 


Limoges & Bergeron (2010) 


12 


56 160 


7.98 


NLTE 


B aimer 


H 


Gianninas et al. (201 1) 


13 


57 000 


7.60 


NLTE 


Lyman 


H 


this work 


14 


55 000 


7.90 


NLTE 


Chandra. 


H+metals stratified 


this work 




-20 -10 10 20 



[A] 

Fig. 15. Fit to the Lyman lines of GD246 (thin line). Overplotted are 
pure H models with r c ff = 57 000 K, log g = 7.6 (thick line), and with 
T cff = 52 000 K, logg = 7.9 (dashed line). 

sets of profiles are rather subtle and the fit quality is similar. We 
verified for some model parameters that the influence of metal- 
line blanketing on the line profiles is negligible. 

We then determined metal abundances. Some of them can be 
compared to the results by other authors (Table 4). 

Carbon No C line is detectable. From the absence of the Cm 
A 1175 A multiplet, we find an upper limit of C/H = 3.2 • 10" 8 . 
All stratified models predict strong C lines, indicating that they 
clearly overpredict the C abundance. 

Oxygen The result for O is similar. We derive an upper limit of 
O/H = 6 • 10 from the absence of a photospheric O vi resonance 
doublet. The stratified models strongly overpredict O. 

Silicon Siiv lines are seen at ^1122.48, 1128.33 A (Fig. 16). 
We find Si/H = 1.2(+0.2) ■ 1(T 7 , that is identical to the result of 
Barstow et al. (2003a). For the diffusion models, the lines fit at 
r eff = 55 000 K and logg = 7.9. 

Phosphorus Fitting the Pv AA 1117.98, 1128.01 A resonance 
doublet (Fig. 16) gives P/H = 4(+2) ■ 10~ 9 . The stratified models 
predict too strong lines. 



Q_ '(/) 03 




!Lt i I i I 

1128.0 1128.4 1128.8 

\[&] 



Fig. 16. Fit to P v and Si iv lines of GD 246 with a homogeneous (thick 
line) and a diffusion model (dashed line). The abundances in the homo- 
geneous model are Si/H = 1.2 • 10~ 7 , P/H = 3 • 10" 9 , and S/H = 5 ■ lO' 9 . 
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Fig. 17. Fit to Ge iv A 1 189.03 A in GD 246 (thin line). Overplotted are a 
homogeneous model (thick line, Ge/H = 5 • \0~~ 9 ) and a stratified model 
(dashed line). 

Sulfur From the absence of the S vi ,U 933.38, 944.52 A reso- 
nance doublet, we find S/H < 5 (±2) • 10~ 9 . The stratified models 
produce excessively strong lines (Fig. 16). 

Iron and Nickel Fe and Ni lines are not detected in the UV 
range. We derive upper limits of 1 ■ 10~ 6 for both species. Fe 
and Ni are overpredicted by the stratified models. 

Germanium Our fit to the Geiv Xll 188.99, 1229.81 A reso- 
nance doublet gives Ge/H = 5(+l) • 10~ 9 (Fig. 17). This is 
about twice the value derived by Vennes et al. (2005) for a 
56 000 K model. The difference can be traced back to NLTE ef- 
fects (Sect. 3.2.1). The diffusion models drastically underpredict 
the line strengths. At = 57 000K and logg = 7.6 no Geiv 
lines are visible. 
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Table 4. Element abundances in LB 1919 and GD246 as determined from UV spectra with homogeneous 
models with 56 000 ±2000 K and log£ = 7.9±0.3forLB 1919 and 57 000 ±2000 K and log g = 7.6 ±0.3 for 
GD 246. Identified ions are noted in the fifth column. Column six indicates the tendency of the abundance 
of a diffusion model to be stronger (>) or weaker (<) than the observed line strength. Literature values for 
the abundances are given in the last column. 



WD 


Element 


Abundance (X/H) 


Uncertainty 


Ions 


Diff. 


Literature 


LB 1919 


C 


< 1.0 


io- 7 


5.0 


io- 8 


Cm 


> 








O 


< 1.0 


io- 6 


2.0 


io- 6 


Ovi 


> 








Si 


< 5.0 


10- 9 


5.0 


io- 9 


Si m+iv 


> 








p 


< 1.0- 


io- 10 


1.5 • 


10 -io 


Pv 


> 








s 


< 1.75 


io- 8 


7.5 


io- 9 


S iv+vi 


> 








Fc 


< 1.0 


10~ 6 


1.0 


10~ 6 


Fe v 


> 








Ni 


< 1.0 


10~ 6 


1.0 


10~ 6 


Niv 


> 






GD246 


C 


< 3.2 


io- 8 


1.5 


io- 8 


Cm 


> 








O 


< 6.0 


io- 8 


5.0 


io- 8 


Ovi 


> 


1.6 


10~ 7 a 




Si 


1.2 


io- 7 


2.0 


io- 8 


Si iv 




5.0 
3.2 
1.2 


10 -8 b 

io- 8 c 

10~ 7 a 




p 


3.0 


10- 9 


2.0 


io- 9 


Pv 


> 


6.3 
7.5 


io- 9 c 

1Q -9b 




s 


5.0 


io- 9 


2.0 


io- 9 


Svi 


> 


< 3.0 


io- 7b 




Gc 


5.0 


io- 9 


1.0 


io- 9 


Geiv 


< 


2.5 


10 -9d 




Fc 


< 1.0 


10~ 6 


1.0 


10~ 6 


Fe v 


> 
> 


< 2.0 

< 1.2 


io- 5b 

IO" 7 a 




Ni 


< 1.0 


10~ 6 


1.0 


IO" 6 


Niv 


> 


< 1.3 


10~ 7 a 



a Barstowetal. (2003b) b Wolff et al. (2001) c Chayer et al. (2001) d Vennes et al. (2005) 
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Fig. 18. Best-fit models to the EUVE spectrum of GD 246 (thin line). 
Top: stratified models; bottom: homogeneous models. Solid lines: mod- 
els including Fe; dashed lines: models include Fe plus Ni. 



type includes Fe plus Ni (setting Fe/H = 1.25 • 1(T 6 and Ni/H = 
7.9 • 10~ 8 in the homogeneous models). Fig. 18 shows that the 
stratified models with Fe as heaviest element provide a better fit 
to the observation than the models also including Ni. A grav- 
ity of \ogg = 8.2 is needed, which is significantly higher than 
the result of the Lyman-line fits. The homogeneous models fail 
completely in reproducing the EUVE SED, whis is in agreement 
with the results of Landenberger-Schuh (2005). 



6.4. Chandra 

Vennes & Dupuis (2002) detected Fe lines in the Chandra spec- 
trum of GD 246. They are identified in Fig. 19. Their NLTE anal- 
ysis indicated Fe/H = 3 • 10~ 7 . 

Fig. 20 displays our best fits using the homogeneous and 
stratified model types as employed for the EfJVE analysis. In 
general, as in the EUVE case, models excluding Ni fit better. 
The final homogeneous model fits worse than the Ni-truncated 
stratified model, and with an unrealistically low gravity of log g 
— 13. The best-fit (Ni-truncated) stratified model gives T e s and 
logg that is in accordance with the most recent NLTE Balmer- 
line analysis of Gianninas et al. (2011); see Fig. 14. From the 
homogeneous models, we confirm the Fe abundance derived by 
Vennes & Dupuis (2002). 



6.3. EUVE 

Two types of models were calculated. Both include the metals 
identified in the FUSE observation (the homogeneous models 
with the determined abundances, and C and O set to their up- 
per limits). For one type of models, Fe is added, and the other 



6.5. Mass and distance 



With r eff = 57 000 ± 2000 K, logg = 7.6 + 0.3, H v = 7.52 ■ 
IO" 4 erg cm" 2 sec" 1 Hz" 1 , m v = 13.09, E(B-V) = 0.0011, we 
derive M = 0.54 + ( J 12 M G and d = 94+^ pc . 
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6.6. Summary on GD246 

Our Lyman-line fit gave T e ff = 57 000 K and \ogg = 7.6. The 
Balmer-line fit by Gianninas et al. (2011) resulted in a simi- 
lar temperature (56 160K) but at a significantly higher gravity 
(7.98). With our parameters, the stratified model compared to 
the FUSE data overpredicts C, O, Fe, Ni, P, and S, but under- 
predicts Ge. That result remains essentially valid for stratified 
models with increased gravity (logg = 7.9). 

No acceptable fit to the EUVE spectrum was achieved. The 
gravity of the best fit homogeneous model is too low (logg = 
7.3), while that of the best-fit stratified model is too high (logg 
= 8.2). It must be stressed that in both cases the formally best 
fits are very poor. 

For fitting the Chandra spectrum, the stratified models are 
only better than the homogeneous models when Ni is arbitrar- 
ily removed. The best-fit stratified model has then log g = 7.9, 
which is closer to the Balmer-line than to the Lyman-line result. 
However, the fits deteriorate when Ni is included. But even when 
Ni is excluded, the best-fit stratified model is not satisfactory. We 
conclude that stratified models do not predict correctly the indi- 
vidual element abundances. This conclusion is corroborated by 
the FUSE analysis. 



7. PG 1520+525 

7.1. Previous investigations 

The first analysis of PG 1520+525 (WD 1520+525) based on op- 
tical spectra was performed by Werner et al. (1991) and gave r e ff 
= 140 000 ± 20 000K and log g = 7.0 ± 1 .0. A first indication that 
the temperature exceeds 140 000 K was derived from the shape 
of an EUVE spectrum (Werner et al. 1996). Utilizing UV spec- 
tra, Dreizler & Heber (1998) found T eff = 150000 + 10000K 
and logg = 7.5 ± 0.5. These parameters, as well as the He, C, O, 
and Ne abundances (Werner et al. 2004b), are given in Table 5. 

7.2. Chandra analysis 

Most counts in the Chandra spectrum are located in the 100- 
120 A region. At longer wavelengths, the ISM absorbs the pho- 
tospheric flux, while at shorter wavelengths, the flux decreases 
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Fig. 20. Best-fit models (thick lines) to the Chandra spectrum of 
GD246. Top: Homogeneous model with abundances according to 
Table 4. Bottom: Stratified model. Both models have no Ni. 
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Fig. 21. Three models with different r cft fitted to the Chandra spectrum 
of PG 1520+525. Model abundances are given in Table 5, logg = 7.5, 
log N(H i) = 20.0. 
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Table 5. Atmospheric parameters of PG 1520+525 and PG 1159-035. 
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Fig. 22. Line identifications in PG 1520+525. Overplotted are two mod- 
els with different T cK : 150 000 K (thick line), 160 000 K (dashed line). 
Model abundances are given in Table 5, logg = 7.5, log N(Hi) = 20.0. 
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Fig. 23. Detail of the Chandra observation of PG 1520+525 (thin line), 
overplotted with two models with solar (dashed line) and ten times solar 
(thick line) Mg abundance (T eff =150 000 K, log g = 7.5). 



due to photospheric O vi bf absorptions. Although the S/N ratio 
is rather poor, it is possible to identify individual lines from O vi 
and Ne vn (Fig. 22). The identification of Mg lines is ambiguous. 
For a detailed list of possible lines in this wavelength range we 
refer to the analysis of H 1504+65 (Werner et al. 2004a). 

We have calculated a small grid of models for T e ft = 
140 000 K, 150 000 K, 160 000 K, and logg = 6.4, 7.0, 7.5. We 
kept fixed the He/C/O/Ne abundance ratio but varied the Mg 
abundance from solar to ten times solar in order to derive an 
upper limit. 

The fit is hardly affected by the particular choice of log g. 
Instead, r e ff is the most sensitive parameter. Mainly judging 
from the spectral slope, we find the best fit at r e ff = 150000 + 
10 000 K adopting logg = 7.5 ± 0.5 (Fig. 22). This confirms the 
previous temperature determination from UV/optical results, but 
the Chandra spectrum allows no further reduction of the error 
bar. 

A change in the Mg abundance from solar to ten times solar 
affects the strength of the lines only marginally, and it is not clear 
which abundance is more likely to fit the observation (Fig. 23). 
The flux of the models is also affected in regions where no Mg 
lines are expected. This can be explained by the occurrence of 
autoionization features in the bf cross sections. These resonances 
can be strong and narrow, looking like a common absorption line 
(Werner et al. 2004a). The wavelengths of these resonances are 
uncertain and could thus be a possible explanation for unidenti- 
fied features in the observed spectrum. 

7.3. Summary on PG 1520+525 

Our analysis of the Chandra spectrum gives an independent con- 
firmation of the result on r ei f derived from previous UV spec- 
troscopy (150000 + 10 000 K, assuming logg = 7.5 + 0.5). 

As outlined in Sect. 2.2, we may use this result to draw con- 
clusions on the location of the blue edge of the GW Vir instabil- 
ity region, which is confined by the non-pulsator PG 1520+525 
and the pulsator PG 1 159-035. The position of both stars in the 
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PG 1520+525: T cff , logg from Dreizler & Heber (1998) and our work. 
Abundances (in mass fractions) from Werner et al. (2004b), except for 
Mg (this work) and Fe (Werner et al. 201 1). PG 1 159-035: from 
Werner et al. (2011). 



r e jf-logg diagram is displayed in Fig. 24. Also shown there is 
the theoretical location of the blue edge predicted by Quirion 
et al. (2007) assuming He/C/O mass fractions of 0.4/0.4/0.2 in 
the stellar envelope, which is within error limits identical to the 
abundances found in both stars (Table 5). Besides these element 
abundances, the exact location of the edge depends on the metal- 
licity; therefore, it is shown here for z = 0.0 and z = 0.007. 
The edge moves to higher effective temperatures with increasing 
z. The z - 0.007 line is appropriate for the two stars consid- 
ered here, because we had found solar Fe abundance in them 
(Werner et al. 201 1). It can be seen that the theoretical position 
of the edge runs between the location of both stars, very close to 
PG 1159-035. 

Also shown in Fig. 24 are other PG 1 159 stars in which pul- 
sations were looked for (taken from the compilation in Werner 
& Herwig 2006). The location of nonpulsators redward of the 
indicated blue edge is possibly due to their different chemical 
composition compared to PG 1159-035, because the composi- 
tion affects the position of the instability strip (Quirion et al. 
2007). 



8. Summary and Conclusions 

We analyzed the Chandra spectra of three hot WDs. Two of them 
are H-rich WDs and the other one is a H-deficient PG 1 159 star. 



8.1. The DA white dwarfs LB 1919 and GD246 

The primary aim of our Chandra observation of LB 1919 was to 
explain its relatively low metallicity compared to similar objects 
as suggested by previous analyses of EUVE spectra. It turned 
out that no metal features were detected in the Chandra spec- 
trum. The same result was found from our analysis of a FUSE 
spectrum. In essence, all data are compatible with the assump- 
tion that LB 1919 has a pure hydrogen atmosphere. This is in 
conflict with our stratified models. In them, the vertical run of 
individual element abundances is computed from the assump- 
tion of equilibrium between gravitational downward pull and ra- 
diative upward acceleration, i.e., the metal abundances are not 
free parameters but computed as functions of T e ff and logg. 
They predict that significant amounts of light and heavy metals 
should be accumulated and readily detectable in the atmosphere 
of LB 1919. 

A few other hot DAs with similar parameters that also do 
not have detected metals are known (Barstow et al. 2003b); 
one famous example is HZ 43 A (r e ff = 51 000 K, logg = 7.9, 
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Beuermann et al. 2006). The reason for the purity of their atmo- 
spheres is completely unknown. It may be speculated that these 
peculiar DAs have no heavy-element reservoir that is assumed to 
be present in the radiative levitation calculations. Barstow et al. 
(2003b) consider whether the possible depletion of those reser- 
voirs by selective mass-loss could be responsible for that phe- 
nomenon or whether the progenitors were metal-poor. They ar- 
gue that, for radiation-driven winds, mass loss should be lower 
in the pure H stars than in those containing heavy elements, be- 
cause the wind is driven by metal lines. Furthermore, these WDs 
are local disc objects, so it seems unlikely that any of the pro- 
genitors could have been metal poor. 

We have also analyzed an archival Chandra observation of 
the hot DA GD 246. Our work was motivated by the fact that 
this is the only DA that exhibits individual metal-line features 
in the soft X-ray range. As for LB 1919, we utilized chemically 
homogeneous as well as stratified models. While we in princi- 
ple expect that the stratified models are a better representation, 
the Chandra observation cannot be fitted satisfactorily by either 
type of model. A better fit by stratified models is only achieved 
when nickel is removed artificially. Otherwise, the atmospheric 
opacity becomes much too large because of a strong overpre- 
diction of nickel. In comparison with the FUSE spectrum, most 
metals are overpredicted, but at the same time, another one (Ge) 
is underpredicted. In the latter case, one can speculate that the 
model atoms are still too small, not having enough line transi- 
tions to absorb photon momentum. The overprediction of most 
metals, however, is more of a problem. We can only speculate 
that an additional physical mechanism, ignored by our models, 
is affecting the equilibrium abundances, e.g., (selective) mass- 
loss or weak magnetic fields. 

The complete failure of our model atmospheres, homoge- 
neous or stratified, to fit the EUVE observation reinforces the 
conclusion that some physics is missing. The same conclusion 
was drawn by Barstow et al. (2003b) from the fact that UV spec- 
tra from DAs with similar T e ff and log g display a wide variety 
of metal abundances. These authors considered the possibility 
that accretion of interstellar or circumstellar matter could cause 
this variety. Recent investigations have revealed the frequent in- 
cidence of dust debris disks around WDs (e.g., Farihi 201 1) that 
are composed of disrupted planetary material. In order to test 
the influence of accretion of such material on the observed pho- 
tospheric abundances, this effect would have to be included in 
the diffusion/levitation models in a manner performed (e.g., by 
Koester 2009) in cooler WDs in which radiative levitation is neg- 
ligible. However, it is difficult to understand how an external 
supply of material would resolve the problem of present diffu- 
sion models generally overpredicting metal abundances. 

Our analysis of the two DAs comprised the derivation of ef- 
fective temperature and gravity from the Lyman lines, and the 
resulting parameters were compared to published results from 
Balmer-line analyses. We confirm previous findings that con- 
flicting results are obtained from the UV and optical line pro- 
file analyses. In particular, Balmer-line temperatures are often 
significantly lower than Lyman-line temperatures (Barstow et al. 
2003a). We found that LB 1919 is another rare case where the 
opposite effect was found. While this obviously points to a short- 
coming in the WD atmosphere models, there is no indication as 
to which physical ingredient is treated inadequately or is miss- 
ing. These uncertainties in the derived atmospheric parameters 
add to the problems in a detailed quantitative comparison, ele- 
ment by element, of observed metal abundances with predictions 
from radiative levitation models. 
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Fig. 24. Pulsating (filled circles) and nonpulsating (empty circles) 
PG 1159 stars in the logr eff - logg diagram. Evolutionary tracks are 
labeled with the respective stellar masses in Mq (Miller Bertolami & 
Althaus 2006). The red edge (short dashed line) of the instability region 
(Quirion et al. 2006) and two blue edges (long dashed lines, Quirion 
et al. 2007) are shown. The upper and lower blue edges are for z = 
and 0.007, respectively. 



8.2. The PG1159 star PG 1520+525 

The Chandra spectrum of the nonpulsator PG 1520+525 was 
used to constrain its effective temperature and to compare it with 
that of the pulsator PG 1159-035. They both confine empiri- 
cally the blue edge of the GW Vir instability region for a par- 
ticular envelope composition. The position of the edge predicted 
by the nonadiabatic pulsation models of Quirion et al. (2007) 
is consistent with the spectroscopic results. This is strong proof 
for the predictive power of these models. The interior structure 
of PG 1159 stars that is inferred by their usage appears to be 
very reliable, so that corresponding asteroseismologic analyses 
are based on solid ground. 
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